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C^- ' Abstract 

o 

r-j ■ A systematic study of the relations between fluctuations of the 

Qh extensive multiparticle variables and integrals of the inclusive multi- 

O 1 particle densities is presented. The generalized factorial moments are 

introduced and their physical meaning discussed. The effects of the 
additive conservation laws are analysed. 



: 1 Introduction 



Studies of multiparticle production is a complicated matter because of the 
large number of variables describing the system and thus several methods 
were developped to simplify the problem. Two of them turned out partic- 
ularly useful in providing information: (i) studies of correlations between 
particles by investigation of the inclusive distributions and (ii) studies of 
event-by-event fluctuations. Since both of them give -in principle- a com- 
plete description of the system, they must be related. This relation was 
firsst discussed in [1] , and reformulated in [2] . In the present paper we follow 
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the arguments of [HE] and investigate further the relation between these two 
ways of approach to multiparticle production. 

To illustrate the problem, consider the well-known formula 



where pk(qi, ••■> Qk) is the inclusive distribution of k particles and F k is the 
factorial moment of the fc-th order: 



F k = Y, P{N)N(N - 1)...{N -k + l) = (N(N - 1)...(N -k + 1)) (2) 



N is the event multiplicity and P{N) denotes the multiplicity distribution. 

The main interest in the formula (pQ) is that its L.H.S. characterizes event- 
by-event fluctuations of the multiplicty while its R.H.S. represents the mea- 
surement of particle distribution by k-arm spectrometer. Thus (pQ) connects 
the two quantities which are defined in an entirely different manner. Needless 
to say, confronting such -apparently unrelated- quantities is often a source of 
a new insight into the problem. This is precisely the interest in investigation 
of further relations of this type. 

Another feature of ([T|) which make this relation very useful in practical 
applications is that it connects the factorial moment of order k with the 
inclusive density of the same order. Thus it can be used even if one does 
not have a complete knowledge of the system (which is of course practically 
always the case)0. 

For k = 2, the extension of ([1]) to quantities other than multiplicity was 
proposed in [3j. It was discussed and applied by several authors (see, e.g., 
[l]-[7] and the review [5]). The possibility of generalization to k greater than 
2 was suggested in [5] and considered in [9]. 

The purpose of the present paper is twofold. First, we derive an explicit 
formula, valid for arbitrary k, which extends the the relation ([1]) to fluctu- 
ations of quantities other than multiplicity. This is obtained by introducing 
the generalized factorial moments [2J, defined in terms of measurements of 
event-by-even fluctuations of measurable extensive quantities. Second, we 
investigate the physical meaning of the generalized factorial moments along 
the lines developped in [10J, as described below. 

1 This feature may be contrasted with the relation between F k and P(N). They must 
be related because they both describe completely the multiplicity distribution. However, 
as seen from ([2]), to calculate F k one needs to know all P(N) (and vice versa). 
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As first discussed in [TU], the factorial moments have a rather simple 
physical interpretation. It follows from the observation that for any distri- 
bution of multiplicity which can be represented as a superposition of Poisson 
distributions 

P(N) = j dNW{N)e~ N — (3) 

one has 

F k = J2( N ( N ~ 1)-(JV — k+ l)P(N) = f dNW(N)N k . (4) 

Thus the factorial moment of the actual multiplicity distribution P(N) mea- 
sures the standard moment of the underlying distribution W(N). One can 
say [10J that the factorial moment removes the statistical noise (represented 
by the Poisson distribution) from the corresponding moment of the distri- 
bution W(N). This observation is easily generalized to local multiplicity 
fluctuations [TO] . 

It is shown in the present paper that, in absence of other constraints, the 
generalized factorial moments can be used as well to remove the statistical 
noise from the data. This is of course to be expected. The problem must 
be treated more carefully, however, in presence of an additive conservation 
law. Indeed, the additive conservation law implies that the "statistical noise" 
cannot be entirely random and therefore the simple argument of [10] must 
be modified. We consider in some detail the effects of conservation of charge 
and of transverse momentum. 

The problem of removing the statistical noise from fluctuations of various 
quantatities observed in multiparticle production, is the long-standing one. 
The most popular method remains the <3> measure, introduced in [TTJ. Other 
possibilities, and their relation to the $ measure were discussed in [5]-[7]. 
They are now commonly used in data analysis [T2]-|17j. 

Our treatment, using the generalized factorial moments, goes beyond the 
previous studies in two points. On the theoretical side, our method applies 
to correlations of any order (until now mostly k = 2 was discussecH). On the 
practical side, including the effects of additive conservation laws makes the 
control over the reliabity of the experimental results much more solicit. 

2 A related method of removing the statistical noise, applicable to arbitrary k, was 
described in [18] and studied in [19] , 

3 For k = 2 the effects of charge conservation were discussed in [6]. An interesting 
proposal to deal with energy-momentum conservation can be found in |20j . 
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In the next section the generalized factorial moments are introduced and 
their relation to integrals of inclusive multiparticle densities explained. Their 
physical interpretation in case of random noise is explained in Section 3. 
Modifications due to conservation of a discrete quantum number are dis- 
cussed on the example of charge distributions in Section 4. The effects of 
conservation of transverse momentum are considered in Section 5. Discus- 
sion and conclusions are given in the last section. Two Appendices explain 
the details of algebra. 



2 Generalized factorial moments 

Generalizations of ([T]) to fluctuations of extensive quantities other than multi- 
plicity were found in [1] and reformulated in [2] . Here we follow the argument 
used in these two papers. 

Consider a single-particle variable x = x(p), where p is the particle mo- 
mentum. Consider, furthermore, the set of extensive quantities 

N 

X l (q 1 ,...,q N ) = Y,Wqn)] 1 (5) 
n=l 

where N is the multiplicity of the event. 

We want to study fluctuations of X = X\. One possible method is to 
consider the moments 

(X k ) = W d qi ...dq N P(q u q N] N) ( Xl + ... + x N f (6) 

N 

where P(qi, g^v; N) is the probability to find an event with iV particles at 
momenta (q%, Qn) and where we have introduced the shorthand 

Xj = x( qj ) (7) 

which will be used henceforth. It was shown in [1] that the moments ([6]) can 
be expressed by linear combinations of the integrals 

R(k lt k s ;s) = J dq 1 ...dq s p s (q 1 ,...,q s )x'l 1 ...x k s s . (8) 

These relations are, however, fairly complicated. 
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To find a more elegant formulation we observe that when one takes x = 1, 
one has X = N and thus the simple relation ([T]) must hold. It is thus clear 
that one has to find a generalization of the factorial moments ([2]). Such a 
generalization was proposed in [2]: 

F k [x] = ([X-(k-l)x][X-(k-2)x]...[X-x]X) (9) 

where x is the operator, acting on a product [X^.-.X^] as follows^ 

i m 

x[x h ...x lm ] = -Y\x h ...x ls+l ...x lm ] (10) 

171 8=1 

where the indices [/i, ...,l m ] need not be different. It follows from (fTUi) that, 
in particular, 

xX l = X l+1 ; x[X,] m = X m [X,] m - 1 (11) 

From (fTU|) one can also deduce how x acts on a product Xj 1 ...Xj k . The 
result is 

1 k 

x[xj 1 ...Xj k ] = {xj l ...Xj k )—'^^Xj a (12) 

K s=l 

where, again, the indices ji, ...,jk need not be different. 
Using © and we obtain for k = 1, 2, 3 

F X [X] = (X) ; F 2 [X] = (X 2 ) - (xX) = (X 2 ) - (X 2 ) ; 
F 3 [X] = (X 3 - 2£X 2 - XxX + 2[x] 2 x) = (x 3 - 3XX 2 + 2X 3 ) . (13) 

We shall now prove that 

N N 

[X-(k-l)x][X-(k-2)x]...[X-x]X = J2 ...J2 x h- X i k ( 14 ) 

»1=1 ik=l 

where all indices are different from each other. 

4 The definition of the operator x given in [2] was incomplete. I would like to thank 
Andrzej Kotanski for pointing out the error. 
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The proof goes by induction. For k = 2 we have X 2 — X 2 = YliLi J2f=i 
with i j. Thus (I14j) is satisfied. Suppose now that it is valid for a given k. 
Multiplying both sides of (fT4l) by (X — kx) we thus have 

N N 

[X - kx][X - (k - l)x]...[X - x]X = (X- kx) ••• J2 x ix- x i k ( 15 ) 

»X=1 *fc = 1 

But 

N N N N N 

X J2 -Y1 x ii- x ik = J2 x il2 ■■■ Z! x ii- x i k = 

h=l i k =l j=l h=l i k =l 

k N N N N N 

= Z J2 ■■■ x^.-Xx^f.-.x^ + Yl ■■■ x u Vis (i6) 

s=lii=l ifc=l s=fe+lii=l *fe=l 

where in the first term in the R.H.S. one of the variables from X is identical 
to one of Xi 1 ...Xi k , whereas in the second term all variables differ from each 
other. Noting that, as seen from f)12p . the first term is identical to 

N N 

kx \Yl ■- Z Zli-ZiJ: (17) 
*l=l *fe =1 

one sees that it cancels in (TH)]) and thus R.H.S. of (TTo^) equals to the second 
term in the R.H.S. of f|T6l) . But this term is just what is needed to complete 
the proof. 

It follows from (TrjJ that 

oo „ 

Fk[X]= E d qi ...dq N P( qi , ...,q N ;N)[X - {k - l)x]...[X - x]X = 

= J p k (qi,---,qk)x(qi)...x(q k )dq 1 ...dq k (18) 

where Pk(qi, Qk) is the fc-particle inclusive density. This follows almost 
directly from definition of inclusive density 

p k { qi , ...,q k )=Y N ( N - 1).--{N — k + 1) ( dq k+1 ...dq N P(q u .., q N ; N) (19) 

N J 

if one notices that the number of terms in the R.H.S. of (1141) equals N(N — 
1)...(N -k + 1). 

Eq. ffl8l) represents a generalization of ([TJ) (which is obtained from fflBl 
by putting x(q) = 1). It has the same two attractive features: it connects the 
fluctuations of the extensive variables Xi with the integral of the inclusive 
density and it involves only moments of finite order. 
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3 Physical interpretation of the generalized 
factorial moments 



To discuss the physical interpretation of the generalized factorial moments, 
we again consider the momentum space split into s-bin and the rest. To 
determine the generalized factorial moment in the s-bin, we need the distri- 
bution of the number of particles and of the variable x in it. Following the 
lines of [10J, we demand that in each bin particles are distributed as ran- 
domly as possible. In absence of any constraints this implies in the s-bin the 
distribution 



77, 

Pa(qi, Qn, n; n, rf)dqi...dq n = e~ n — J[ [g(qi; rj)dqi] (20) 



n[ i=i 



where n is the average multiplicity in the s-bin and g(q; rj) is the momentum 
distribution of one particle on this bin, with r\ representing a collection of 
parameters on which this distribution may depend. 

To obtain the actual distribution we have to weight (j20~|) by the underlying 
probability distribution W s (fi, rj): 

P s [qi, q n , n] = J dndrj m W s (n, rf)p s (qi, q n , n; n, rj). (21) 

From (12TT) we deduce the inclusive particle density in the s-bin 

Pk(qi, ■~qk) = ^2n s (n s - l)...(n s - k + 1) / dq k+1 ...dq n P s (qx, ...,g„;n) (22) 

n s J 

Introducing (|22|) into (|18|) and using fTSU]) we obtain for the generalized fac- 
torial moment in the s-bin 

F k [X] = J2 n ( n - - k + 1) / dq 1 ...dq n P s (q 1 , ...,q n ;n)x(q 1 )...x(q k ) 

n J 

= J dndr]W s {n,r])[nx] k = ([nx] k ) = (x fc )(23) 

where 

x = x{rj) = J g(q; rf)x{q)dq (24) 
is the average of x(q) at fixed rj. X is the average of X at fixed n and r\. 
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Eq. (I2"5j) is the generalization we were seeking for. It shows the physical 
interpretation of the generalized factorial moments: they remove the statis- 
tical noise from the moments of fix = X of the underlying distribution W s . 
For x = 1 we of course recover the well-known result [10J. 



4 Fluctuations of charge 

As already explained in the introduction, the argument present in the previ- 
ous section fails in case when the variable x satisfies an additive conservation 
law. In this section we discuss fluctuations of charge. The discussion applies 
to any discrete, additive quantum number. 

We again select an s-bin from the momentum space and consider the 
distribution of particles inside and outside of this bin. 

The first problem is to define the distribution describing the "statitical 
noise" . Following the idea of [10] , we take the normalized product of Poisson 
distributions (expressing the independent particle emission and thus intro- 
ducing no correlations). This is supplemented by the Kronecker ^-symbol to 
satisfy the conservation law: 



Ps = i-e — -r-e — — -e -— e 



A Q n+! n~\ N*+\ N*~\ 

8n++N*+-n--N*--Q (25) 

where n + , n~ denote the number of produced charges in the s-bin, iV* + , iV* - 
those outside of the s-bin. Q is the net charge in the initial state. We take 
Qo > 0. The normalization factor Aq is given by (see Appendix 1) 

A Q = e-*u>- Q °I Qo {N) (26) 
where Iq denotes the Bessel function, 

A> = fi + + N* + + fT + N*- =N + + N- (27) 

and 

rz — ~ 2^" 
N = 2^N+N- ; uo = = (28) 

N 2N+ 
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The observed distribution of particles is thus 



P(„+,„-,N~,N-;Q ) = J < B- < B- d N-UN-W^,n^,N-) 

p.(n + , n ;N'+,N-;Q ; n+ , ft" ; iV* + , N'~) (29) 

where is the "dynamical" distribution, free of the statistical noise^]. 

It is important to observe that, contrary to the naive expectation, the 
parameters n + ,h~; N* + , N*~ do not represent the average values of n + ,n~; 
N* + ,N*~, respectively. The explicit formulae for these average values are 
given in Appendix 1. In particular, the average values of the number of 
produced positive and negative particles in the the s-bin are given by 



i W(#) \ = / n± N I QoTl (N) \ 
I Qo (N) I \N±2 I Qq (N) I 



(30) 



where we have denoted by~the average at fixed h + ,n~; N* + , N*~ , and by (...) 

the average over n + ,n~; N* + , N*~ (with the probability distribution W). 
The explicit formulae for the generalized factorial moments of the charge 
and the (standard) factorial moments of the multiplicity in the s-bin 

are derived in Appendix 1. They read: 




fc! 



' (k — m)\m\ 



n 



N+ 



k—m 



n 



+2m-k 



(AO 



Iq (N) 



(31) 



Eqs. (1301) and fl3TT) simplify substantially in the interesting limit of very 
large number of produced particles, N ^ N — > oo. In this limit, using 
the asymptotic expansion of the two Bessel functions, one obtains up to the 
terms of order N~ l 



F ± 
r k 



\n+ ± n-] k - k{k ^- l) [n + ± n-f- 2 



2N 



(32) 



The first term corresponds to the standard interpretation of the general- 
ized factorial moments, as expressed by ( 1231) . One sees from ( )32i) that the 



The selection of statistical noise in the form of (|25jl is not unique. Another natural 
possibility is to consider independent emission of particle pairs (supplemented by indepen- 
dent emission of Qq positive particles). The two methods differ mostly by the treatment of 
the initial charge Qq. Without further information about dynamics it is difficult to decide 
which of them describes better the physics of the problem. 
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correction to this result (induced by the conservation law) vanishes, at fixed 
n , with the inverse power of the total multiplicity of produced particles. In 
this case the correction is expected to be small at high energies (particularly 
for the central heavy ion collisions^. We thus conclude that the generalized 
factorial moments can indeed provide a useful tool for eliminating the sta- 
tistical noise from the event-by-event fluctuations of multiplicity (Fj^) and 
charge (Fj^) even in presence of the conservation law. 



Fluctuations in transverse momentum 



In this section we discuss fluctuations of the tranverse momentum, as an 
example of a continuous variable subject to an additive conservation law. 

Selecting one bin in rapidity (s-bin) and dividing the available momentum 
phase-space into this bin and the rest, we write the transverse momentum 
distribution in the form 

P[q u Qn, n; q* 1} q* N *,N*; q ] = 
dndqdDdN*dq*dD*W(h, q, d, ; N*, q*, D*) 



n 



III 



n 

1=1 



-N' 



N*\ 



N* N* 

n 

1=1 



2tt(D 



*\2 



— %l + ... + q n + q{ + ... + q* N * 
A i 



(33) 



where q denotes the two-dimensional transverse momentum vector and go is 
the total transverse momentum of the considered phase-space region. The 
uncorrelated " statistical noise" in the form of Poisson distribution for multi- 
plicity multiplied by a gaussian distribution of transverse momenta is supple- 
mented by the 8 function insuring the transverse momentum conservation^ 
and by the normalization factor A±. 
It is also natural to take 



hq + N*q* 



(34) 



6 For F k the relative correction may be large if n+ + n_ 3> n+ 



7 This is a commonly used form of the "statistical noise" (see, e.g., [2Ql [21] ). Note 
that, since W(n,q,d,;N*,q*,D*) is a general (positive) function, the assumed gaussian 
forms in (|33|) do not restrict seriously the observed distribution of particles (as both the 
average values, q, q* and dispersions, D, D* of the Gaussians can fluctuate according to 
the distribution W). 
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Replacing the S function in (1331) by its Fourier transform we obtain for 
the inclusive distribution in the s-bin 

p k ( qi ,...q k ) = J dNdqdDdN*dq*dD*W(n, q, D; N*, q* , D*) 

eN ^ q y-n y /2 _ 1} 1 dye ^y q o^n(n-l)...(n-k + l)e 



(2tt) 2 A ± J » V n! 



dq {k+1) ...dq n J] 



1 3 igij/„-ta--?)72D 2 



(35) 



Using ([TBI we thus have for the generalized factorial moment in the bin s 

Ski. _ I ' s k 



F fc [X;g ] = / dndqdDdN*dq*dD*W(n,q,D;N*,q*,D*)^ = (-^) (36) 

J On I \ On I / 



with 

,fe± = (2tt) 2 
(note that v4j_ = Sq±), and 



= [ dye*™ [n<P(y)f e^ m ~ D '^'-%^'- D * V/2 -H (37) 



0(2/) = ^5 / rf 2 ge-^) 2 / 2D2 e^a;(g) ; 0(0) = x (38) 

where 5 is the average value of x at fixed q and D 2 . The function </>(?/) 
depends obviously on the choice of the variable x(q). Generally, if x(q) is a 
polynomial in q of order I, then (f)(y) is a polynomial in y of the same ordeJ§ 
multiplied by e *?W^72. 

If the number of particles is large, the integral over d 2 y can be evaluated 
by the saddle point method. We write 

S k± = -^Jd 2 y[n<f>(y)] k e^ (39) 

with 

= -i 9oy + fi[e i(iy - D2y2/2 - 1] + jV*[ e ^-°*V/2 _ y ( 40 ) 



8 E.g., for x(g) = <? we have 0(y) — (q + iyD 2 )e^ v y D I 2 , for x(q) = q 2 we obtain 
4>{y) = + (q + iyD 2 ) 2 }e^y-y 2D ' 2 / 2 . 
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The saddle point equation \1>' = gives y ~ q — nq — (N*q*) and thus 
(IMj) implies y = 0. Furthermore \I>"(0) = -N{D 2 + q 2 ) = -NA 2 , where 
JV = h + iV* and 

~ 2 nD 2 + N*D* 2 _ 2 n<f + iV*g* 2 

-D = ^ ; q = — — — (41) 

N N 



Consequently we obtain 
S ki NA 2 



Sq± 2n 



d 2 y[n^y)] K e-y 2NA2 / 2 . (42) 



To evaluate this integral in the limit of large N, one can expand 4>(y) in 
powers of y and observe that, since the coefficient in the exponent increases 
with increasing N, the dominant term at large N is that with the lowest power 
of y 2 under the integral (the terms with odd powers of y do not contribute). 
A detailed discussion is given in Appendix 2. Here we anly summarize the 
results. 

If 0(0) = i^0 one obtains 



F k [X) = {[nq)") (43) 

i.e. we recover the formula (1231) . Thus, if the average value of x(q) in the 
s-bin does not vanish, the standard interpretation of the generalized facto- 
rial moments remains valid, even in presence of the conservation law. The 
corrections to this result vanish as 1/ < N >. They are given in (|62|) of 
Appendix 2. 

This is not the case if x = and qx ^ 0. For k even, k = 2p, one obtains 

F 2p [X] = (-lfp\ {[nqxY (44) 

which shows that now the factorial moments are related to the average value 
qx rather than to x. For a fixed n and large N these moments tend to zero. 
For a finite ratio n/N, however, they may be large. The corrections to ([41 
and the formula for k odd are given in Appendix 2. 



6 Discussion 



We have reconsidered relations between the event-by-event fluctuations of 
extensive multiparticle variables and the integrals of the inclusive distribu- 
tions [1J. Our main result is the formula 018p which expresses the integral 
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of the fc-particle inclusive distribution in terms of the generalized factorial 
moment of the k-th order. 

The physical meaning of the generalized factorial moments is discussed. 
It is shown that they remove the random uncorrelated statistical noise from 
the data, the result known already from previous investigations [TQl [2] . When 
an additive conservation law is at work, however, the statistical noise cannot 
be uncorrelated and, consequently, its removal can at best be approximate. 
The exact formulae were derived for charge and transverse momentum con- 
servation. The corrections were evaluated in the limit of very large total 
number of produced particles, relevant for collisions at high incident energy. 

In short, the results presented in this paper formulate a systematic ap- 
proach to investigation of fluctuations of extensive variables. They seem to 
be particularly useful for studies aiming to uncover the local structure of the 
multiparticle system. To take just one example, interpretation of the data on 
transverse momenta (or energies) in small rapidity bins will be much more 
transparent if presented in terms of the generalized factorial moments. 

It was shown in []Jj that the relation between fluctuations and correlations 
can be extended to other extensive variables. This, however, goes beyond the 
scope of the present investigation. 



7 Appendix 1 

Replacing the Kronecker 5 symbol by its integral representation one obtains 
P(n + ,n-;iV* + ,;V*-;Q ) = (w Q {n + , n~; N* + , N*~; Q )) (45) 

where 

1 1 



WQ (n + } n-;N* + ,N*-;Q ) = — — [ dfe 

An 2ir Jo 



iQof 



n+\ n~\ N*+\ N*~\ { } 

and Aq is the normalization factor. 

Thus the generating function for the distribution wq is 

<^ Q (z + ,z-,Z + ,Z-) = 

E w Q (n\n-;N*\N*-;Q )(z + T + (z-r(Z + ) N+ *(Z-) N -* = 

n±,iV±* 

= J_ df e -iQof e n + (z+e>f-l) e n-(z-e->f-l) e N+*(Z+e*f-l) e N-*(Z-e-'f~l) 
27T JO 
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From this it is not difficult to derive, by the standard methods, the 
multiplicities. One obtains 

1 r27r 

fj± = _JL_ f dfe -iQof e N+(e if -l) e N-{e-«-l) h ± e ±if = 



2nA Q Jo 

2mA Q J z 

00 ~ 1 r dz 
e -—n ±l u ±l -^ Y UN)— l— z iz ±l ~^ 



where oj and N are defined in 
Similarly, one obtains 



It follows that 

, N N ra ± 

The inclusive distribution in the bin s is 

where 

H Q (k + ,k-,f) = e** + ^ f - 1 '>e ff *~^ f - 1 '> V(n + ,n-;k + ,k-) 

rt+ ,n~ 



rr r \ n~\ 

and 

k\ n + \ n~\ 



V(n + ,n-;k + ,k-) = 5, 



k~k+~k- 



k+\k-\ (n+ - k+)\ (n~ - k~)\ 



— dhp~ ihk k\ H+] r ihk+ H ' r ihk ~ 

2nJo k+\(n+-k+)\ k-\(n--k-)\ 

14 



is the number of ways one can select k + out of n + positive and k~ out of n~ 
negative particles in a given order. 

To calculate the generalized factorial moment in the s-bin 

Fk= E (-l) k ' Pk (k + ,k-) (54) 

k++k~=k 

we observe that 

V Y j (-l) k ~V(n + ,n-;k + ,k-) = — / dhe~ ihk \\ + e ih } n [1 - e ih ] n (55) 
Tt u 2tc Jo 



k+ k 

and thus 



3 7V+(e^-l) p 7V-(e-^-l)_^_ /" 27r ^ £ -ihk £ [h+ -h~ e -^]e ih 



27T 

= e ^ + (ev-i) e ^-(e-v-i)^+ e i/ _ ^- e -i/]fc ( 56 ) 

Now integration over df: 

1 f 2jT 



fATT 



k++k~=k 

— [ 2n dfe-We^'-^-^'-Vin+eV -h~e^} k 



2ir Jo 



6 f (IjZ AT 



-e 



N(z+l/z)/2 { 



ujz) Qo [h + ujz-n /{ujz)} h 



2ni J z 

k k\ 

N,.-Qo £ \-h- /uj] m [n + uj} k - m 



e u) 



m =o (k-m)\m\ 

1 r dz 



j=-oo 

fc! 



iV '-" Q ° E 71 ,y my [-n-M m [un+} k - m I Qo+2m ^ k (N), (57) 



e uj 



m=0 



(A; — m)!m! 



where A r± , iV and are defined in Eqs. fl27|) and (|28|) . 

The normalization factor Aq is obtained from (jSTj) by putting k — 0. 
Inserting ([57} and ([50]) into (|5IJ) and (J53D we obtain f[3Tj) . 

We also note that the standard factorial moment for the multiplicity can 
be obtained from f[5Tj) simply by replacing [— h~] m by [ra - ]" 1 . 
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8 Appendix 2 



Here we evaluate the generalized factorial moments for variables related to 
transverse momentum, given by (1421) . in the limit of very large total number 
of produced particles (but without any restriction on the ratio n/N, i.e. on 
the size of the s-bin). 

The starting point is the formula . 

Developping the function <f)(y) around y = we have 

m = E tm^r (58) 

where 

K = / d 2 qx(q)q m e-^/ 2D2 = 4^ (59) 

This implies that [4>{y)] k can also be represented by a series 

oo 

[<P(y)] k = E € } y m (so) 

m=0 

where the coefficients 0w can be explicitely evaluated, when necessary. 
Introducing this in (PEZl) we have 




One sees that in the limit of very large N the first nonvanishing term domi- 
nates. If 4>q = x we have 

F kX [x(q)) « (n*#) = ([nx) k ) (62) 

and thus we recover the original formula (1231 . We conclude that in this 
case the standard interpretation of the generalized factorial moments remains 
valid even in presence of the conservation law. 
If x = 0, and 0i = qx ^ 0, we have 

= [iy<p 1 } k + (iy) k+l <p k l - 1 <p 2 /2 + .... (63) 



[<Kv)Y 



[iyfa] k 



20i 
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For k even, k = 2p the dominant term is 



r^IW (fgY (64) 



where we have used the relation ( l59l) . 
For odd, fc = 2p — 1 we obtain 



= (-If hp-\qxr 2 ^p\ {j^J ■ (65) 

where again flBTJl) was used. Note that these moments are proportional to q 2 x 
and thus vanish if the variable x(q) is odd with respect to change (q <-> —5). 

We conclude that the standard interpretation of the factorial moments, 
as expressed in (1231) . holds only when the average value x does not vanish. 
When x = the factorial moments measure the average value qx rather than 
x. 
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